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TECHNICAL NOTE 2794 


A COMPARISON OF TWO METHODS OF LINEARIZED 
CHARACTERISTICS FOR A SIMPLE 
UNSTEADY FLOW 
By Roger D. Sullivan 


SUMMARY 


Two methods of using the concept of linearized characteristics are 
derived for the one-dimensional unsteady flow in a tube that is rotated 
about an axis perpendicular to the axis of the tube. One of the methods 
corresponds to that used by Ferri in his basic work on the subject. 
Solutions are made .by both methods for boundary conditions that allow 
analytic solutions. Comparison shows that both methods give the same 
results but there are significant differences in their application. ^ 


INTRODUCTION 


The term linearized characteristics has been applied by Ferri in 
reference 1 to a process of superposing a small perturbation on a basic 
solution of a set of nonlinear hyperbolic differential equations. The 
perturbation can be due to a change in the prescribed conditions along 
the boundary, in the position of the boundary, or in the differential 
equations themselves. The unknowns can each be expanded, for example, 
into a power series in a parameter representative of the perturbation, 
with unknown coefficients, except for the coefficient of the first term, 
representing the basic solution, which is known. The equations for the 
other coefficients are found to be linear and can be solved by charac- 
teristic methods. As applied to steady supersonic flow, this procedure 
has the advantage, stated roughly, over ordinary perturbation theory 
that it can be used for perturbations of any known flow (to which the 
characteristic method of solution can be applied) rather than only for 
perturbations of constant parallel flow. 

A version of the process was used by Sauer (ref. 2) for determining 
the steady supersonic flow about bodies of revolution at small angles of 
attack but apparently has not been used for other problems except in 
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reference 1, where. its application to a variety of _problems is discussed. 
In references 1 and 2, the expansion of the unknowns was substituted 
into the original equations of motion, which were then transformed into 
characteristic form. The purpose of the present paper is to compare 
that-method of procedure (called method A subsequently) with an alter- 
nate method (method B) whereby expansions of the dependent and inde- 
pendent variables are substituted into the characteristic form of the 
equations. , ........ 

The comparison is made by developing the equations of both methods 
as they apply to a simple set of nonlinear hyperbolic equations of gas 
dynamics, the set relating to the one-dimensional unsteady motion of~a 
gas in a tube of constant cross section with the addition of a perturba- 
tion due to the rotation of the tube. The equations derived for this 
case are then solved for a set of boundary conditions for which the 
solution can be found analytically by both methods. 


SYMBOLS 


a, A 
c v 

f,g 

k 

P 

R 

s ,S ■ 

t, T 

u, U 
x,X 
a,p 
7 

e 


speed of sound 

specific heat at constant volume 

unknown functions of one variable 

constant 

pressure 

gas constant 

entropy increment divided by gas constant 
time 

velocity 

distance along tube measured from axis of-rotation 
characteristic parameters 
ratio of-specific heats 
absolute- temperature 
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R 7-1 

P density 

cd angular velocity 

Subscripts: 

Letter subscripts indicate partial derivatives. 

Number subscripts indicate order of approximation with the following 
exceptions: 

Sq reference value of a 

xq and xj_ inner and outer ends of tube, respectively 

ANALYSIS OF UNSTEADY ROTATING FLOW 
Development of Equations 


The linearized-characteristics concept is applied to a relatively 
simple problem In order to bring out more clearly some of its essential 
features. Consider the unsteady motion of an inviscid gas in a cylin- 
drical tube which can“rotate~ with angular velocity to about an axis 
perpendicular to the generators of the cylinder (fig. 1). Distance 
along the tube Is denoted by x measured from the axis of rotation. 
Only values of x large in comparison with the greatest lateral dimen- 
sion of the tube in a plane of rotation are considered so that the 
local velocity of the rotation may be taken as cdx. Cross flows and 
lateral pressure gradients are neglected so that the flow considered 
Is one-dimensional. A perfect gas with constant specific heats is 
assumed. 

When CD = 0, the ordinary equations of motion apply (ref. 3); that 
is, the continuity equation 

p t + PUx + up x = 0 (!) 

and the momentum equation 

u t + uu x = " £ Px 


( 2 ) 
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When o> ^ 0, equation (l) is unaltered, hut equation (2) has another 
term representing the acceleration along the tube due to the rotation 


u t 


+ UUv 


2 

ox 



(3) 


The entropy is now assumed tro be constant - throughout the gas so that 


p = kp 7 


(4) 


See appendix A for a treatment of the case of variable entropy.' Equa- 
tion (4), the equation of state 


p = pRS 


and the relation 


a 2 = = 7R9 

dp 


allow the replacement of p and p by functions of a in equations (1) 
and (3). Thus the equations of motion become 


where 


la-(. +- aUjj- + Xua x = 0 


u^ + uu x + Xaajj - cu 2 x = 0 


(5) 


X = 


2c v = 2 : 

R 7-1 


These equations are referred to subsequently as the "primary equations" 
to distinguish them from the characteristic equations which are now 
developed. 


Adding and subtracting equations (5) gives 


u^ + (u + aJujj +— xja^ + (u + a)a^J - o) 2 x = 0 
u^ + (u - a)u x - X ja^ + (u - a)a0 - o^x = 0 


( 6 ) 


By the usual procedure for hyperbolic equations (ref. 3)> new independent 
variables a and 3 are introduced so that one equation of equations (6) 
has derivatives with respect to a only and the other, with respect 
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to (3 only. Then the lines a = constant and f3 = constant form the 
two sets of characteristics. This process is accomplished by writing 


*a = (u + a)t a 
xp = (u - a) tp 


( 7 ) 


so that 


= t a 
oa, 


a + 


+ (u + a) ^ 


St 


Sx 


A + (u - a)A 

St Sx 


Thus, equations (6) become 


Uq, + Xa a - a>^xt a = 0 

up - Xap - o^xtp = 0 


( 8 ) 


Equations ( 7 ) and (8) form a set of four equations in the four unknowns 
u, a, x, and t and are equivalent to equations (5)* These equations 
can always be solved by step-by-step methods. 


developed. Let 


Meth od A . - Linearized approximations to equations ( 5 ) are now 

u(x,t) = U + CD^U-^ + + • ■ . 

a(x, t) = A + cr^a-L + cs^ag + . . . 


(9) 


Only even powers of co are used since equations (5) show that the problem 
is independent of the sign of o>. By substituting expressions (9) into 
equations (5) and setting the coefficients of each power of cd equal to 
zero, the following equations are obtained: 


XA t + AU X + AUA X = 0 


U t + UU X + XAA x = 0 


(10) 
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Xai-fc + Au]_ x + XUai x + U x a]_ + XA x ui =- 0 
u lt + ^ Uu lx + XAa ix +-- U X U]_ + XA x ax - x = 0 


( 11 ) 


Xa, 


2t 


+ Au£ X + MJa2 X + U x a2 + X^i^ + a l u lx + *- u l a lx = 0 


u 2t + U u 2x +" ^ a 2x ‘ K ^x u 2 + lA-x a 2 + u l u lx + ^l a lx = 0 


( 12 ) 


and so forth. Or, when these equations are put into, characteristic form 
by the same method that was applied to equations (5), the same equations 
for the characteristics are found to apply to each setj that is. 


Xo = (U + A)t a 
xp = (U - A) tp 
whereas the other equations are successively 

u a + XAq, = 0 
dip - XAp 


( 13 ) 


- 


( 14 ) 


u la + ^la + + ^x)ta( u l + a l) “ x ^a - 0 

u lg - Is-ip + (U x " ^x)^p( u l ~ a l) - x tp = 0 

+ Xa 2a + ( Ux + M x) t a( u 2 + a 2) + ( u lx + Xa lx) t a( u l + a l) = 0 

Ugp - Xagp +- (Ux _ ^X ) "k(3 ( U 2 a 2) + ( U lx ~ l a lx)^p( u l _ a l) = ® 

and so forth. 


( 15 ) 


(16) 


The development and use of equations ( 13 ) to (16) is referred to 
hereafter as "method A." 


Method B.- Instead of"' substituting expansions of u and a into 
the -primary equations (eqs. ( 5 )) as in method A, these expansions are 
now substituted into tine characteristic form of the equations (eqs. (8)). 
But the values u and a must be considered functions of a and 3 
in equations (8), and the corresponding values of x and t are found 
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from equations (7). Hence, x and t are also expanded. The method 
is, therefore, similar to the technique presented by Lighthill in refer- 
ences 4 and 5* 


By writing 

u(a,3) = U + cd^IJ^ + cd *»2 ♦ • • • 
a(a,3) = A + cd^Aq_ + cn^Ag + . . . 


( 17 a) 


x(a,3) = X +■ co^X}. + 01 ^X 2 + • • • 
t(a, 3 ) = T + cd%2. + ^^2 + • • • 


(1713) 


substituting into equations ( 7 ) and (8), and setting the coefficient of 
each power of cd equal to zero, the following sets of equations are 
obtained: 


U<x + — ® 


1 


Uo 


AAp — 0 


Xo - (U + A)T a = 0 


X p - (U - A)T p = 0 


u la + ^loc 
u lf3 “ ^13 



(18a) 


(18b) 


(19a) 


X ia - (U + A)T la = (Ui + A x )T a 
X lp - (U - A)T lp = (U-L - A l} Tp 


( 19 b) 




U 2a + ^2a = X l T a 

U23 - XA2JJ = X-jTp 


+ XT 


la 


+ XT 


13 


(20a) 
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X 2a - (U + A)T 2a = (Ui + Ap)Tp a + (U 2 + A 2 )T a ] 
X 2 p - (U - A)T 2p = (U X - Ap)Tpp + ( U 2 -^ A 2 ) T p| 


(20b) 


and. so forth. The development and use of these equations is referred 
to as "method B." 


Equations (l8) and equations (13) and ( 14 ) are-exactly the same 
and their solution is the complete solution of the problem when co = 0. 
This solution, the basic solution, is essential to the solution of the 
higher approximations and is considered' known. The equations for the 
higher approximations in both methods are all linear in the unknowns. 


Application to Specific Boundary Conditions 

In order to demonstrate some of the differences between the two 
methods, an application to a simple case for which the equations can 
be integrated analytically is now made. Assume that the tube extends 
from x = x 0 to x = Xj > x 0 ) and is closed at the ends. The 

tube is initially at rest, and the state of the gas is given by u = 0 
and a = a Q , a constant. At time tr = 0 the tube is started rotating 

about x - 0 with angular velocity co. The neglect of cross flowB and 
lateral pressure gradients might be completely unjustified in this case; 
nevertheless," the problem is interesting mathematically. 

The problem, then, is to solve the primary equations (eqs. ( 5 )) 
with the following boundary conditions: 

At- _ t = 0 , 

u = 0 ; a = Sq (Xq ^ x < xp) ( 21 a) 


At x = xq and x =■ xp, 


u = 0 


(t >0) 


The basic solution (i.e., for co 


0 ) is, of-ncourse , 


U = 
A = 



(21b) 


( 22 ) 
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Since these are constants, this problem is similar to that of ordinary 
linearized theory for steady supersonic flow about slender bodies. 


The introduction of a and' 3 is somewhat arbitrary as long as 
equations ( 13 ) are satisfied; these parameters may be conveniently 
taken so that the two families of characteristics of the basic solution 
are given by 


x = a^t + 3 I 
x = -aQt + a 


(23) 


so that 


x = g-tl 
2 



> 


( 24 ) 


For method B, x and t are replaced by X and T in equations (23) 
and (24) . 


When cd has a finite value, the xt-plane may be divided into 
regions by the characteristics which pass through the points (xq,0) 
and (X]_,0) as shown in figure 2. The solution in region I is 

unaffected by the walls at the ends. In region II the effect of the 
wall at x = xq only is felt; whereas, in region III that at x = x^ 

only is felt. Both walls affect the flow in region IV. 

The solution in region I is found first. (in appendix B the 
solution in region I is found in closed form. ) With method A, equa- 
tions (22) and (24) are substituted into equations (15)- The first 
equation of equations ( 15 ) is integrated with respect to a, the second 
with respect to 3 > so that 

u i + 

"1 - - - 5= 


a 2 + 2a3 + f( 3 ) 


2a3 + 3 ^ + g(c 


( 25 ) 
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where f((3) and g(a) are arbitrary functions subject - to the- -boundary 
conditions. These functions are determined by noting that at t = 0, 
a, = p, by equations' (23) , and that U]_ - a^ = 0 by equations ( 21 a). 

It is found that f(p) = -33 2 and g(a) = -3a 2 . By substituting these 
expressions into equations ( 25 ) and solving for uj_ and a^, the fol- 
lowing expressions are obtained: 



( 26 ) 


Equations (l 6 ) can now be solved in the same manner. Thus, 


U 2 = - 


3a 0 3 




> 


(27) 


®2 = 


2X +- 3 fa - p\^ _ ( 2X + 3)a 0 ^ 


24\ 2 a 0 3 


Substituting these expressions into equations (9) gives as the solution 
in region I ' 


2 1 k o 

U = 03 Xt - i 03 Xt j + 


— = 1 - — 03 2 t^ +- — — t . 3 . 03^^+ . . 

a 0 2X 2kx 2 _ 




( 28 ) 


The solution by method B starts with equations ( 19 a), which are 
the same as equations (15) in this case. The solution, then, is 



(29a) 
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Now that these values are known, equations (194) may he solved. The 
arbitrary functions are determined by making lines in this solution 
pass through the same points on the x-axis of the xt-plane as the cor- 
responding lines of the basic solutienj that is, X]_ = Tq = 0 when 
a, = p . Thus , 




> 


(29b) 


= 4 + 1/tt - P\ 3 _ x + 1 t 3 
"6Xa Q 3V 2 j ‘ 61 


Solving equations (20) in the same way gives 


U 2 = 2k. 


IX + 1/q + p\ (a - B y = 2X. + 1 

6X a 3V 2 A 2 ) ex 


A2 — 


■0 

2.X 




24A 2 


=-^-^a 0 T 
S-r\3 ' 2 ' 24X 2 


(30a) 


X 2 


T 


2 


+ 4 / q, + p\ (a - = 5A + 4 4 

24^ 2 A 2 y " 24*.^ 

(X + 1)(5A + 7)/a - p\ 5 = (x + lHjjX + 7 ) t 5 
120A 2 a 0 5 ' 2 / 120X 2 


( 30 b) 


Substitution into the series (17) gives as the solution for region I 


u 


cjq 2 XT 


+ 


Si-JU, aha? 

6x 


+ • • 


_S_ = 1 - _i_ co2t 2 - g A + 1 ^4^4 + 
a 0 2X 2U 2 


n 




(31a) 


u 
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x = X + i a&CT 2 + cd^CT^ 4- -. . . 

2 2k\ 

t - T + — - cd 2 ! 3 4 - ^5 + 

6 \ 120 Xr 


> (31b) 


The variables X and T may be interpreted simply as parameters 
of the solution. Actually, they define points of intersection of* -char- 
acteristics in the basic solution corresponding to characteristics which 
intersect at x, t * in. this solution. In any case, X and T may be 
eliminated from equations ( 31 ); the procedure results in equations ( 28 ) 
again. 

The only difference in the solution for region II is in the deter- 
mination of the arbitrary functions resulting from integration. Those 
functions resulting from integration with respect to p may be deter- 
mined by the condition that the solution must ‘equal that for region I 
on the boundary between the solutions j that Is, when p = Xq. The other 

functions can then be found from the boundary condition given by equa- 
tions (21b). That condition, however, does not apply to the solution 
of- equations” (19b) and (20b). For these equations, the boundary con- 
dition is replaced simply by 


x = X +- 


co 2 X n 


+ CD 




. = Xq 


or 


X - XqJ 


x x = X 2 = . . 


. = 0 
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The solution found by method A is 


u = ||( x - x o)( a o t + x o) - gI ^ ~ (x “ x o)^ 8Xa o 3t3 + 


3(3^ - 7)xQa 0 2 t 2 - 12j7>>. + k)xQ - (X + 2)xjxoaQt 

Rl7X + 3l)x 0 - (5X + 7)xj (x - x Q ): 


,x or + 


55- = 1 - ^f~2&v t2 + 2 W - 2x o( x - *oj 


CD 


24x 2 a Q 4 


(2\ + 3)a Q 4 t 4 + (3X + 5)x Q a3t3 + I2(x - 


3 J^( A. + 3)xq - (5X •)- 7)xj (x - x 0 )x 0 a 0 t - 
4(X 4- 2)(4xq - X)(X - Xq) 2 Xq| 4- . . . 


>-(32) 

x o) x o a o 2t2 + 
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and by method B 


u = SL 


Xa 0 T - XQ(a 0 T - X + xq) 


cs 


4— 


-1 24Xa Q 3 


x^aoT - X + Xo)j4(2\ + l)ao 2 T 2 - (13X + 11) (X I xq^T 

( x - xo)|u - i)x + (nx + i3)x^j-j + . . . 


4(2X + l)xa 3 t3 


1 


>(33a) 


_a _ or 


*0 


2Xa r 


a 0 ^ 2 + 2xo(a 0 T - X + xq)] ^-j-((2X +-l)a Q V 

— 24X^a r 4 1 


x 0( a 0 T - x + xq)<(7X + 


5)a Q 2 T 2 - + (] 


(11X +-13)x n - (5X + 7)X 


IX a Q T 


2(X - xq) 


(5X + 7)x n - 2(X + 2)X 


X 


M + . 


..2 I - 


x = X + 


0) 


2Xa,_ 


2 m 2 


XXa Q T - x 0 (a 0 T - X + x 0 ) ( Xa Q T 4- X - jcq) 


t T + ^ R* + ’^ x Q^ a 0 T "■ x + x o)( a O T +■ x ~ * 0 ) 

12Xa 3 - - -J 


> (33b) 


4 * . • 


The solutions for region III are obtained in a '.similar way and are 
the same except that Xq is replaced by xj_ and the signs of some terms 

are reversed. In regions II and III, it is found, that the elimination 
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of X and T from the solution by method B results again in the solu- 
tion by method A. 

Figure 3 shows characteristic diagrams for the case = 2xq, 

03 = i ^2-, and X = 5 ( X = lA). The dashed lines are the character - 

2 Xq 

istics of the basic solution and, therefore, are the mathematical 
characteristics used in the higher-order solutions; these lines are in 
the xt-plane for method A and in the XT -plane for method B. The solid 
lines are the true characteristic lines of the problem in the xt-plane. 
They are found automatically by method B. Method A, as described so 
far, does not give the true characteristic lines, but it is found that 
the solution of equations ( 7 ).> after the substitution of the values of 
u and a (for example, eqs. ( 28 ) for region i), results in the same 
characteristic lines as have been found by method B. This determination 
is necessary for the evaluation of the flow in the field, as illustrated 
in figure 4 where typical distributions of velocity and speed of sound 
are shown for a given time and in figure 5 where typical variations with 
time at a given value of x are shown. In figures 4(a) and 5(a), the 
results found by method A without the determination of the variation of 
the characteristics are shown. That is, the lines PR and PS of 
figure 3 are used as the boundaries between the regions. It is seen 
that u and a are discontinuous at these boundaries (although the 
quantity u - Xa is continuous along PR and u + Xa is continuous 
along PS). In figures 4(b) and 5(b), the results are shown when the 
true characteristics of the problem, the lineB QR and QS of figure 3* 
are taken as the boundaries. Method B automatically gives the results 
of figures 4(b) and 5(b). Parts (a) and (b) are identical except where 
they represent conditions in the regions between PR and QR and 
between PS and QS in figure 3« 


DISCUSSION 


Some aspects of the preceding example which are due to the fact 
that U and A are constants and which would not ordinarily be present 
should be noted. (l) The unsolved equations of method A are unusually 
simple in comparison with those of method B since the terms containing 
the unknowns explicitly drop out. (2) As a result, the solution of 
the equations of’ the first approximation is the same by both methods. 

( 3) The solution in method B for Xq_ and T-j_ affects the unknowns 

u and a only in the second order and similarly for the higher 
approximations . 

The basic solution cannot usually be found analytically but must be 
determined by a step-by-step procedure. In that case, the approximations 
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must also be done step by step. The characteristic equations would be 
written in a finite difference form with Ax and A3 cancelled out. 
Thus, a and p would not be used explicitly. With method B, the 
calculations would be carried out in the XT -plane where the character- 
istic directions remain constant as they do in the xt^-plane in method A. 
The transfer to the xt-plane, which consists of merely adding the terms 
of equations ( 17 b)., can be delayed until the desired number of approxi- 
mations has been made. 

From the preceding simple example, some inferences can be drawn 
concerning the application of the two methods to other basic equations. 
The unknowns do not occur explicitly (although they do appear as deriv- 
atives) in equations ( 19 a) and (20a), but, if they occur thus in the 
basic characteristic equations, corresponding to equations (l8a), they 
will occur in the equations corresponding to (I9a) and (20a). If the 
original independent variables (X and T in the given example) appear 
in the basic characteristic equations, the four equations in each approx- 
imation of method B must all be solved simultaneously, rather than as 
two pairs of equations as in' the example. This complication is a serious 
disadvantage of- method B and may rule it out completely in such cases. 

The relative simplicity of the equations of method B, taken one by 
one, is probably quite general. The appearance of derivatives of lower 
approximations with respect to the original independent variables in the 
equations of method A is to be expected as long as the primary equations 
are nonlinear, whereas they cannot appear in method B. 

If an application is made to the two-dimensional supersonic flow 
about a single body, for example, and conditions on the surface of- the 
body only are desired, the calculation of the true characteristics is 
superfluous. In cases such as this, method A has the advantage of having 
two less equations to solve in each approximation. However, the pre- 
ceding problem indicates that the true position of characteristics 
through points of discontinuity in the boundary conditions is necessary 
for the determination of the flow in the field and for the determination 
of conditions on another boundary if such characteristics intersect it. 

This necessity for correcting characteristics as it^applies to 
linearized theory for steady supersonic flow is discussed In reference 6. 
There, as in the present example, what is here called the basic solution 
consists of constants so that this effect is felt only in the second and 
higher approximations. In more general cases it would be felt in the 
first approximation. 

In problems where the perturbation involves a change in the position 
of the boundary, method B has the advantage that the position is constant 
in the plane of the original independent variables of the basic equations 
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(for example, the XT -plane in the present problem) where the solution 
is carried out. 

When the perturbation involves the addition of a new dependent 
variable, it can be handled in much the same way as the entropy is 
handled in this case (appendix A). 


CONCLUDING REMARKS 


Two methods of using the concept of linearized characteristics are 
derived for the one-dimensional unsteady flow in a rotated tube. A com- 
parison of both methods has shown that they give the same results but 
there are significant differences in their application. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 11, 1952. 
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'APPENDIX A 

DEVELOPMENT OF EQUATIONS FOR VARIABLE ENTROPY 

When the entropy is not constant throughout the fluid, equations (5) 
are replaced by 


la-t + au x + Xua x = 0 


u-t + uu x + \aa x - a 2 s x - o> 2 x = 


where s is the entropy increment divided by the gas constant. As long 
as the flow is without shocks, the entropy obeys the relation 

st + us x =0 (35) 

The characteristic equations become 



and equation (35) may be replaced by 


s a tp + spta - 0 (37) 

(See ref. 3 for a derivation of the preceding equations for cu = 0.) 
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Writing 


s(x,t) = S + oi^S]_ + 0i^S2 + . . . 

substituting this expression along with equations (9) into equations (3*0 
and. (35), and converting to characteristic form give the equations of 
method A as 


Ucc, + XA<x - — ASa = 0 


U 3 - XA p + i AS p = 0 > 


S to + S R t_ = 0 
a p pa 


( 38 ) 


•la + 

Xaia + ( 

; Ux + 


a i) 

■IP " 

Xa lp + ( 

; u x - < 

1 

1 — 1 

ea 

-P 

a i) 

la'tp 

+ s ip t a 

+ 2u 1 S x t a t p = 0 




( 39 ) 


+ Xa 2a + ( U x + ^x)t a (u 2 + a 2 ) + (u lx + Xa lx )t a (u 1 + a J 




J^2a + (2a X + u l) s la + ( 2a 2 + «2) s a + ( U 1 + a l) 2s x t ^j = 0 
u 2p - Xagp + (U x - XaJ t p (112 - a 2 ) + (u lx - ^i x ) t p( u l ' a l) 
^[as 23 + (2a x - u a ) s ip + (2^ - u 2 )s p - (u x - a^) 2 S x t p J = 0 



(b0) 


s 2a t P + * 2 ^ 0 . + 2 ( u l s lx + u 2 S x) ^p = 0 
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and so forth. Equations (13) still give the. mathematical characteristic 
direction for each approximation. If the true characteristics are to be 
found, equations (7) still apply. 

In order to obtain the equations for method B, the expression 
s(a,3) = S +- a3^S]_ + oj^S 2 + • • • 

along with equations (17) is substituted into equations ( 36 ) and (37)> 
so that ■ -- - 


U a + XA<x - y AS a - 0 

Up - *Ap + y ASp = 0 

Xa - (U +-A)T a = 0 
Xp - (U - A)Tp = 0 
Sal 1 3 + SpT a = 0 


“A 


(hi) 


u la + 

^■la 

- ~(AlS a 

+ . AS la) 

= xt a 

%3 - 

\A3_p 

+ ^( a 1 s 3 

+ AS ip ) 

= XTp 


Xia - (U +■ A) T la = (Ui +■ Ai)T a 
X ip - (U - A)T 13 = (% - A x )Tp 

SiaTp + SaTip + SipT a + SpTi a = 0 


(te) 
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U 2a + XA 2a “ y( A '2 S a + A l S la + AS 2a) “ X l T a + ^ la 
U2(3 - *A 2 |3 + A-2 S P + A 1 S 1P + AS 2p) = X 1 T P + 
x 2a " (U + A)T 2a = (U X + AjT^ + (u 2 + A 2 )T a 
X 2 p - (U - A)T 2 p = (Ux - Ax)t 1P + (U 2 - A 2 )t P 
s 2a.Tp + S la Tip + S a T 2 p + S2pT a + S^Tia + SpTpa = 0 


and. so forth. 

It is to he noted that, when the entropy of the primary solution S 
is a constant, the equations of both methods simplify considerably. For 
each approximation in this case, the equations for the entropy may be 
solved first, then the equations for the velocity and speed of sound, 
and then (in method B) the characteristic direction equations. Other- 
wise, all three equations (in method A) or all five equations (in 
method B) in each approximation must be solved simultaneously. 
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APPENDIX B 

SOLUTION IN CLOSED FORM FOR REGION I 


The form of equations (28) suggests the- substitutions 

u = xf(t) 
a = a 0 g(t) 

Equations (5) become 

Xgf + gf = 0 (44) 

+ - f 2 - co 2 = 0 ' (45) 

at 

whereas the boundary condition given by equations (21a) becomes 


f(0) = 0; g(0) =1 

Equation (45) is easily integrated to give 

f = co tanh cat 

By substituting this value for f, equation (44) can be solved tor-give 

g = ( cosh cot) 

Therefore ' . - — 

u = cox tanh cot 

— = (cosh cot) 

SO 

These functions may be expanded in power series; the procedure results 
in equations (28) once more. 
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Figure 3-~ Characteristic lines for the case 

X = 5 (7 « 




X 1 = 2 * 0 ' “ = | ^ 8114 
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(a) With mathematical characteristics of method A as boundaries between 

regions. 



1-0 1.1 1.8 1.3 1A 1.5 1.6 1.7 l.fi 1.9 8.0 

JL 

*0 


(b) With true characteristics as boundaries between regions. 

i ant 

Figure 4 .- Distribution of velocity and speed of sound when — — =0.3 

x 0 

T ®“0 

for the case Xn = 2 xq, q = -x — , and X = 5 (7 = 1 A) . 






